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J% ; Abstract 

Finite temperature Euclidean SU(2) lattice gauge fields generated in the con- 
(— i \ finement phase close to the deconfinement phase transition are subjected to cool- 

ing. The aim is to identify long-living, almost-classical local excitations which 
carry (generically non-integer) topological charge. Two kinds of spatial bound- 
! ary conditions (fixed holonomy and standard periodic boundary conditions) are 

applied. For the lowest-action almost-classical configurations we find that their 
relative probability semi-quantitatively agrees for both types of boundary condi- 
tions. We find calorons with unit topological charge as well as (anti-) self dual lumps 
(BPS-monopoles or dyons) combined in pairs of non-integer (equal or opposite sign) 
topological charge. For calorons and separated pairs of equal-sign dyons obtained by 
cooling we have found that (i) the gluon field is well-described by Kraan-van Baal 
solutions of the Euclidean Yang-Mills field equations and (ii) the lowest Wilson- 
fermion modes are well-described by analytic solutions of the corresponding Dirac 
equation. For metastable configurations found at higher action, the multi-center 
structure can be interpreted in terms of dyons and antidyons, using the gluonic and 
fermionic indicators as in the dyon-pair case. Additionally, the Abelian monopole 
structure and field strength correlators between the centers are useful to analyse 
the configurations in terms of dyonic constituents. We argue that a semi-classical 
approximation of the non-zero temperature path integral should be built on super- 
positions of solutions with non-trivial holonomy. 
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1 Introduction 



It is widely accepted that quark confinement in QCD is related to some complex struc- 
ture of the gauge field vacuum. A popular scenario views the vacuum state as a dual 
superconductor |]J in four dimensions which is unable to sustain strong gluon electric 
fields such that flux tubes are formed due to a dual Meissner effect. The analogon to the 
Cooper pairs of usual superconductors are Abelian monopoles, in other words, condensed 
magnetic currents belonging to an Abelian subgroup of the strong gauge group. 

The lattice evidence in support of this mechanism rests either on the ability to localize 
the magnetic currents 0, |3| as the worldlines of magnetic defects with respect to some 
gauge fixing prescription M, or on the construction of a corresponding monopole creation 
operator || in order to study its condensation. The question how their condensation 
leads to confinement is relatively well-understood || [7J, and while also the question of 
universality with respect to the gauge condition is partly answered affirmatively || U , the 
reason why these monopoles are created in the QCD vacuum is far from being understood. 



There is another working picture of the Yang-Mills vacuum JTOl |TT|, O, |13| which 



has been applied very successfully to hadron physics, including basic features like chiral 
symmetry breaking and the ?7a(1) anomaly as well as details of spectroscopy []14[]. It 
is based on the instanton, a solution [15| of the Euclidean field equations. However, 
even in the form of the instanton liquid model [IB], when some local interaction between 



instantons is taken into account, instantons cannot be brought into relation with the third 
basic feature of QCD, confinement. There have been attempts to stretch the instanton 
liquid model to its limits in order to provide a confining string tension of sufficient strength 
between a static quark and antiquark ]|TTfl , at least at intermediate distances of a few 
instanton radii. A recent assessment of the instanton-generated forces between static 
charges has been discussed in Ref. 



In lack of a satisfactory instanton-based mechanism of confinement, many indications 
have been presented that instantons (or, more generally, carriers of topological charge) 
are closely related to the Abelian monopoles which are usually detected by gauge fixing 
procedures. This evidence stems from the observation made for one-instanton and many- 
instanton configurations [IS, |2D| in the maximally Abelian gauge and from the observation 



of short range correlations between topological charge and magnetic currents in genuine 
confining lattice configurations |2l], |22| . On the other hand, also in artificial instanton liq- 
uid model ensembles (multi(anti)instanton configurations) monopole network percolation 
has been observed p3f as a prerequisite of monopole condensation. This has led to the 
conclusion that semi-classical, smooth gauge field configurations can give rise to networks 
of (light, condensed) magnetic monopoles even if these become discovered only as defects 
of the gauge fixing process. Thus, some reconciliation of the two pictures seems possible. 

However, as long as the instantons are uncorrelated, the Abelian monopoles and cen- 
ter vortices [Q revealed by the respective gauge fixing are quantitatively insufficient to 
provide confinement [|ll|. Moreover, the approximate Casimir scaling at intermediate 
distance that the strong force acting on charges of different representations obeys is also 
strongly broken in the instanton liquid ensemble ||18||. This makes it less likely that in- 
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stantons are related to confinement even if some important correlation would be included 
into the model. 

Nowadays, in a fresh attempt to relate instanton physics to the confinement property 
of the Yang-Mills vacuum, instantons dissociated into meron pairs have been investigated 
with respect to their monopole and center vortex content p5| . Our present work, which 
also aims to extend the instanton picture, starts from a somewhat different corner, finite 
temperature. All the criticism with respect to an instanton mechanism of confinement 
applies also to finite temperature Yang-Mills theory below the deconfinement temperature, 
where periodic instantons (calorons) are the ingredients of the instanton liquid model 



More precisely, there are periodic instantons with trivial holonomy which are thought of 
being the bricks of topological structure. An isolated caloron with trivial holonomy is a 
periodic classical solution of the Euclidean equations of motion with trivial asymptotics 
of the Polyakov loop at infinity, P(x — > oo) = z G Z(N). These are the "old" calorons 
first considered by Harrington and Shepard [E7|. 

During the last years new caloron solutions have been found and studied by Kraan 



and van Baal (KvB) |29[ which correspond to a non-trivial asymptotic holonomy 
P(x — > oo) G" Z(N). From our point of view, a particularly interesting feature is that, in 
the generic case, these configurations are composites of BPS-monopoles 0] (or dyons) 
i. e. lumps carrying non-integer topological charge. A semi-classical description of the 
(finite temperature) vacuum in the confinement phase, based on these new calorons still 



remains to be developed. Some new perspectives have been pointed out in Ref. pi| . In 



the present paper we investigate metastable configurations which appear in the process of 
cooling [B2| in the region of few times the instanton action. These are excitations which 
might become the building blocks of such a semi-classical description. The topological 
content, which appears at this deep, almost classical level of cooling, is itself characteristic 
for the phase (confining or deconfined) where the configurations are taken from before they 
are cooled. Still, we do not claim that the number and size of topological lumps obtained 
in this way immediately characterizes the finite-temperature vacuum. 

The outcome of cooling with normal periodic boundary conditions has been compared 
with that of cooling with spatial boundary conditions of prescribed holonomy, which is 
enforced by cold boundary links appropriately chosen. As for the confinement phase, 
the dependence on the boundary conditions is weak. A closer check shows that, after 
cooling with periodic boundary conditions, the average Polyakov loop over regions of low 
action and topological density is not restricted to trivial holonomy. In other words, in the 
confined phase non-trivial holonomy boundary conditions for topological charge lumps are 
provided by the dynamics itself. They do not need to be enforced. In the confined phase 
these are the natural fluctuating asymptotic conditions for the semi-classical excitations. 

Similarly to the studies of the monopole and vortex content of meron [25] pairs, we 
study here also the monopole content and other signatures of the simplest of these con- 
figurations (to be interpreted as dyon-dyon and dyon-antidyon pairs) as potential semi- 
classical building blocks of a confining vacuum. Preliminary results have been presented 



in 031, 34, 35 



In contrast to the meron studies fl25 |, instead of tayloring particular semi-classical 
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configurations, in the present work we are going to demonstrate the emergence under 
cooling of metastable lattice configurations with low action which (distinguished by action 
plateaus) resemble the new calorons (dyon-dyon pairs). In this case, where an analytical 
solution is available, we show that action density, topological density and the profile of 
the fermionic zero modes on the lattice can be fit simultaneously with the corresponding 
expressions for the new caloron KvB solution |28|, [36] . It is worth noticing that at actions 
around one instanton action also configurations pop up under cooling which resemble 
dyon-antidyon pairs. This interpretation is suggested by action and topological density 
and is corroborated by the zero modes of the Wilson fermion matrix and the Polyakov loop 
profile. Lacking an analytical solution for these configurations of mixed (anti) self duality 
it is impossible to present a fit of the various distributions. 

Also for higher metastable plateaus with an action of few times the instanton action a 
multi-center structure of few-lump dyon-antidyon mixed configurations can be recognized 
having a broad distribution of holonomy away from the lumps of action and topological 
charge. We do not claim that this directly proves that the finite temperature vacuum is 
composed exactly as that gas of dyons and antidyons. Nevertheless, it is remarkable that 
cooling starting from genuine Monte Carlo configurations in the deconfined phase ends in 
completely different low-action configurations as it has been already claimed in an early 
related work by Laursen and Schierholz |37 |. 



The paper is organized as follows. 

In Section ^| we will collect some necessary formulae covering the zero-temperature 
instanton, the well-known self dual finite temperature instanton (the "old caloron"), the 
new KvB calorons, monopoles as instanton chains etc. This will be needed for the fit of 
dyon-dyon pairs becoming visible on the lattice. 

Section || provides all necessary lattice definitions, in particular the boundary condi- 
tions and the observables considered in order to identify KvB and other solutions. 

In Section ^ we describe the results of a first part of our simulations where equilibrium 
configurations are cooled down to an action near to the one-instanton action. Although 
we cool to this low level of action, we find the topological content (and the dependence on 
the boundary conditions) very different, depending on whether the original Monte Carlo 
ensemble was describing the confinement or deconfinement phase. 

In Section |5| we report on a second part of our study, where higher lying action plateaus, 
appearing in the process of cooling with periodic boundary conditions, are identified with 
respect to their dyon content. 

Section || presents the conclusions. 

2 Calorons with trivial and non-trivial holonomy in 
continuum SU(2) Yang-Mills theory 

Throughout this paper SU (2) gauge theory in four-dimensional Euclidean space is consid- 
ered. We start the description of caloron solutions with generically non-trivial holonomy 



with the well-known one-instanton solution WB, In the singular gauge its gauge 
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potentials are given by 

A,(x) = ^r ] ; u d u log<j ) (x), 0(a) = 1+ P " v = 1, • • • , 4 (1) 

where 77" is the 't Hooft tensor, r a /2 - denote the generators of the SU(2) group, p, xq are 
the size of the instanton and the 4D position of its center. For further use we introduce 
the following notations 



(p{x) 



ip(x) 

tP{x) = 2rr 2 ((x-x ) 2 + p 2 ), ij>(x) = 2n 2 (x - x ) 2 , (2) 



and rewrite the instanton fields in the form 



A li (x)= -r]^T 3 d u log (j)(x) 

+ \<P{^)Re ((77^ - ifjl^in + iT 2 )d v x{x)) . (3) 

The next step is to construct the time-periodic instanton or caloron [|27|]. For this purpose 
we generalize the function in Eq. (JJ) 

n=oo 2 

m = 1+ E 



((f-x ) 2 + (t-nb) 2 ) 

rrp 2 sinh(gf) 
&r cosh(^) - cos( 



(4) 



where r = \x — x^\ and t = x<± — x\ . b is the time period of the periodic instanton 
solution. For simplicity, b = 1 is assumed in the following. Physically, the periodicity b 
has to be identified with the inverse temperature T _1 . Then the caloron potentials are 
given by Eq. (|3|) with 

<p{x) = , 

ip{x) 

TVO 2 

ipix) = cosh(27rr) — cos(27rt) H sinh(27rr) , ip(x) = cosh(27rr) — cos(27rt) , (5) 

r 

1 rep 2 sinh(27rr) 
x (x) = 1 - - = - 



This solution can be viewed as a chain of instantons aligned along the temporal direction, 
separated by distance b but with the same orientation in color space. In the limit p —>■ 00 
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it can be transformed into a static solution by an appropriate gauge transformation. This 



static object can be identified with a BPS monopole |30[ in Euclidean space, where the 
fourth component A4 plays the role of the Higgs field in the adjoint representation. The 
solution has both electric and magnetic charge. Therefore, we call it also dyon (D). 

Constructing a caloron with nontrivial holonomy p8f is the next step. It can be 



approximately viewed as an instanton chain with periodicity b, too, where each of the 
instantons is rotated with respect to the previous one by an angle 4ttuj in color space, 
uj being the parameter of holonomy as we will see below. The rotation axis can be any, 
for definiteness let us take the third one. The words approximately viewed mean that 
when the instantons are well separated (p/b « 1) the fields near the instanton centers 
look approximately as described above. The caloron field with nontrivial holonomy is 
described again by Eq. (^) but now with pg | 



</>{x) 



ip(x) 
tp(x) 



( r 2 + g 2 +7r 2 4\ 

ip(x) = — cos(27rt) + cosh(47rra;) cosh(47rso;) H sinh(47rra;) sinh(47rsa;) 

2rs 

+7rp 2 (s -1 sinh(47rsu;) cosh(47rrk)) + r _1 sinh(47rra;) cosh(47rsco>)) , (6) 

(r 2 + s 2 -7c 2 p 4 ) 

ip(x) = — cos(27rt) + cosh(47rrco>) cosh(47rsu;) H sinh(47rra;) sinh(47rsu;) , 

2 

x (x) = e A7Tituj ^- {e^V 1 sinh(47rscj) + r" 1 sinh(47rrcu)} 

instead of Eq. ([5]). The holonomy parameters uj and uj are related to each other as 
u = 1/2 — uj, < uj < 1/2. r = \x — x\\ and s = \x — X2I are the 3D distances to the 
locations of the two centers of the new caloron solution, 7ip 2 /b = d = \xi — x 2 \ is expressed 
by the distance between the centers. From Eqs. (^|) it can be easily seen that, when 
uj — » or uj —>■ 0, the caloron with nontrivial holonomy turns into the Harrington- Shepard 
caloron described by Eqs. (^). When the separation between the centers becomes large, 
d = np 2 /b » 1, two well-separated constituents emerge which are static in time. The 
"mass" ratio of these dissociated constituents is equal to uj/uj. Since the full solution is 
selfdual, the ratio is the same for the action as for the equal-sign topological charge carried 
by the constituents, the latter summing up to one unit of topological charge, Q t = 1. The 
separated constituents form a pair of BPS monopoles (or dyons) with opposite magnetic 
charges. In the following we will call it a DD pair, while non-dissociated calorons will be 
denoted as CAL, although all these objects represent limiting cases of one and the same 
solution. The single dyon originally obtained from the Harrington- Shepard caloron in the 
infinite size limit p — ► 00 can be recovered from the new solution by sending the mass of 
the second constituent to zero and simultaneously its position to spatial infinity. 

The action density in all three cases described above can be expressed by a simple 
formula 

^) = -7^>g^(x) . (7) 
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where 



So far, the new caloron solution is in the so-called algebraic gauge . It can be made 
periodic by a gauge transformation which is non-periodic in time g(x) = e - 2mtuJ - T 'i } 

\nl v rzd v \og(f> (8) 

+ 2 ^ Re K^i" ~ if llu)( T i + ir 2 )(d u + 4niuj5 uA )x) + <W 27ru;r3 , 

2 

~ = e -4^ x = r e -2mt s -i sinh ( 47r SLJ ) + r -i S inh(47rrcj) } . (9) 

Now the time component of the caloron potential has become nonzero at spatial infinity. 
We can define the holonomy which becomes non-trivial 

V{x) = P exp (i A 4 (x,t)dt \ -> = e 2wiuJT3 for \x\ -»• oo. (10) 

In terms of cj, the normalized trace of the holonomy, the Polyakov loop which we shall 
take as a direct measure of the holonomy, at spatial infinity becomes 

L(x) = -tr V(x) -> Loo = ^tr 7^ = cos(2vtcj) . (11) 

In the special case = u; = 1/4 the measure of holonomy is equal to zero and the two 
constituent dyons acquire equal mass, i. e. equal action and topological charge. For our 
later purposes it is useful to remember that the Polyakov loop has peaks L(x) = ±1 very 
close to the center positions x = x\^ of the constituents |28| . 



For well-separated dyons, when the functions 4>(x) and ip(x) are almost time indepen- 
dent, the strongest time dependence comes from the first part of the function x(x). This 
dependence is represented by the phase e~ 2mt and is nothing else but the homogeneous 
rotation of the first dyon, which has L(x[) = —1, in color space around the third axis with 
angle 2n over the period. The second dyon with L(x 2 ) = +1 is static. Such a relative ro- 
tation of two dyons (that form a monopole-antimonopole pair) gives the so-called Taubes 
winding necessary two produce unit topological charge from a monopole-antimonopole 
pair d39[. One can detect this rotation in a gauge invariant fashion by looking at the 
gauge invariant field strength correlator defined on each constant-time slice and watching 
its evolution over the periodicity interval h. The field strength is selfdual (E% = B^. or 
antiselfdual, E% = —BD everywhere in the KvB caloron. The electric field E%(xi, t) at the 
centers of both dyons % — 1, 2 is proportional to an orthogonal matrix (in both the SU(2) 
color and space indices). Thus the three (k = 1,2,3) components E%(aii) of the electric 
field form vectors in color space which represent a local reference frame at the center of 
the first dyon. The comparison with the local frame at the center of the other dyon, x 2 > 
can be made in a gauge invariant manner by connecting the centers by the (fixed time) 
Schwinger line parallel transporter 

S(x 1} t]x 2 ,t) = P exp (^i j A£(aT',t)y dx' k ^ . (12) 
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Using this one can form the gauge invariant field strength product 



R X ki{t) = tr ( E k (xx,t) S(x h t;x 2 ,t) Ei(x 2 ,t) S + {x 1 , t; x 2 , t) ) 



(13) 



which is again an orthogonal matrix. This matrix performs a full rotation with the 
Euclidean time t running from to b. 

Finally, let us comment on the zero-mode eigenfunctions of the fermionic massless 
Dirac operator in the background of the KvB solutions. They have been studied analyti- 
cally in f4(| and f36j. One finds closed solutions depending on the type of (anti)periodic 
boundary conditions (b.c.) imposed on the fermion fields in the imaginary time direction. 
In case of well-separated dyon pairs, i. e. for d = up 2 /b » 1, the zero eigenmode 
densities become very simple expressions 



|Vr(x)| 2 

|^ + (x)| 2 



— ^—d 2 [tanh(27rrt<;)/r 

4:71 



4.7T 



d 2 [tanh(2™) /s) 



for antiperiodic b.c. 
for periodic b.c. . 



(14) 



This means that the zero-mode eigenfunctions are localized always around one of the 
constituents of the KvB solution, for antiperiodic b. c. at that constituent which has 
L(x) = — 1 at its center X\. Switching to periodic b.c. for the fermion fields the zero 
mode localization jumps to the other constituent monopole of the gauge field. Therefore, 
the fermionic zero modes provide a convenient way to identify a monopole-pair structure 
in the gauge fields. 



3 Detecting dyons and calorons on the lattice 

Our first aim was to detect the simplest dyon configurations in the context of finite 
temperature lattice simulations. For this purpose we have considered SU(2) lattice gauge 
theory on an asymmetric lattice using the standard Wilson plaquette action with coupling 

S = ^2s(x,t) = ^2^s(x, t;/J,,v), (15) 

X,t X,t fJ,<U 

s(x, t; n,v) = p (1 - ~tr U Xiia> ), U Xtflv = U Xilt U x+ftiV Ul +f>il pl iV 

and periodic boundary conditions in all four (toroidal) directions. For simplicity the 
lattice spacing is set to a = 1. The lattice size will be x N t with the spatial extension 
N s = 16 or 24 and with T _1 = Nt = 4. For N t = 4. the model is known to undergo the 



deconfinement phase transition at the critical coupling [3 C ~ 2.299 |gT|. Throughout this 
paper we concentrate on the confinement phase, i.e. j3 < (3 C . 

We shall generate the quantum gauge field ensemble {U XtlJ ,} by simulating the canonical 
partition function using the standard heat bath Monte Carlo method. The equilibrium 
field configurations will be cooled by iteratively minimizing the action S. Usually, cooling 
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in one or another form is used in order to smooth out short-range fluctuations, while 
(initially) leaving some long-range physics intact. The cooling method applied here is the 
standard relaxation method described long time ago in p2| . 

This method, if applied without any further limitation, easily finds approximate so- 
lutions of the lattice field equations as shoulders (plateaus) of action in the relaxation 
history. Under certain circumstances, this defines and preserves the total topological 
charge of a configuration. However, the short-range structure of the vacuum fields is 
changed. Still, the type of classical solutions, which are being selected, depends on the 
phase which the quantum ensemble {U Xjtl } is meant to describe. We want to investigate 
smoothed fields at different stages of cooling, by using a stopping criterium which selects 
the plateaus in a given interval of action. First, in Sect. 4, we focus at the lowest ac- 
tion plateaus, i.e. S = (0.5- •■ 1.5) S inst , where S inst = Sit 2 / g 2 . denotes the action of a 
continuum instanton. Lateron, in Sect. 5, we shall describe more complex approximate 
solutions found by stopping at higher plateaus. 

The smoothed fields will be analyzed according to the spatial distributions of the 
following observables: 

• action density computed from the local plaquette values and averaged with respect 
to the time variable: 

* t 

• topological density computed with the standard discretization and averaged over the 
time variable: 

Qt(x) = -^ 2 4 -32tt 2 ^ ( S ^pati[U x ^ u Ux, P a]) ; (17) 

• Polyakov loop defined as: 

1 Nt 

L{x) = -ti\{U x - Ml (18) 
t=i 

where the Us 7 t,4 represent the links in time direction; 

• Abelian magnetic fluxes and monopole charges defined within the maximally Abelian 
gauge (MAG). The latter is found by maximizing the gauge functional 

FM = ^tr(^r 3 ^ r3 ) ; (19) 

x,n 

with respect to gauge transformations U x ^ — > C/| = g(x)U x>IJi g* (x + ft). Abelian 
link angles 9 Xtfi are then defined by Abelian projection onto the diagonal U (1) part of 
the link variables U x>fl G SU{2). According to the DeGrand-Toussaint prescription 
[0 a gauge invariant magnetic flux P through an oriented plaquette p = (x, \xv) 
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is defined by splitting the plaquette Q p = 9 x>fl + 9 X +^ V — 9 x+ p^ — 9 XiV into 6 P = 
Q p + 2Tin p , n p = 0, ±1, ±2 such that 9 P G (— 7r, +tt}. The magnetic charge of an 
elementary 3-cube c is then m c = ^ ^ pG a c @p ! 

• eigenvalues and eigenmode densities of the non-Hermitean standard Wilson-Dirac 
operator 

2J D[U) xra>ys p i) s p{y) = A ^ ra (x) (20) 

with 

studied both with time-antiperiodic and time-periodic boundary conditions. For 
our purposes it will be sufficient to consider this operator which breaks explicitely 
chiral invariance. To use a chirally improved lattice Dirac operator would be a next 
step. We find the A spectrum and the eigenf unctions with the help of the implicitely 
restarted Arnoldi method and use the standard ARPACK code package for 

this aim. 

For production and cooling of the equilibrium gauge field configurations we shall use two 
kinds of spatial boundary conditions as in |34j and with preliminary results presented in 

1. standard periodic boundary conditions (p.b.c) on the 4D torus; 

2. fixed holonomy boundary conditions (f.h.b.c): fixed holonomy is realized by cold 
timelike links Ug jt> 4 on the spatial boundary Q. 

For clarity we stress, that the second case is periodic, too, but for the spatial boundary 

tt = {x I x = (l,x 2 ,x 3 ), (xi, 1,2:3) or (x 1 ,x 2 , 1)} 

all time-like links Ugf,i are frozen to constant group elements. For definiteness we have 
used embedded pure Abelian link variables = cos 9 + ia 3 sin 9. In the confinement 

phase at (3 < (3 C we require L(x) = = (L) (corresponding to holonomy parameter 
uo = 1/4) which is satified by 9 = ir/2N t . As in we have studied also the deconfinement 
case ((3 > (3 C ). In this case we fixed the boundary time-like links such that L(x) = (L) 
for x G £1. In both cases (L) denotes the ensemble average < | L(x)\/Ng > of the 
volume- averaged Polyakov loop. 

Each kind of boundary conditions will be employed both for the Monte Carlo produc- 
tion of configurations and for their subsequent cooling. 



D[U) 



xra,ysfi 
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4 Dyonic lumps and other objects seen on the lattice 



In a first part of our investigation we have searched for topologically non-trivial objects 
with lowest possible action, late in the cooling history, in order to find systematic depen- 
dences of the selected solutions on the spatial boundary conditions and on the tempera- 
ture of the original Monte Carlo ensemble. Cooling was stopped at some (n-th) cooling 
iteration step when the following criteria for the action S n were simultaneously fulfilled: 

• S n < 1.5 Sinst, 

• \S n — S n -i\ < 0.01S'i ns t, 

• S n — 2 o n _i + S n -2 < 0. 

The last condition means that cooling just passed a point of inflection. 

For each (3 < (3 C we have scanned the topological content of 0(200) configurations. In 
this late stage of cooling we find approximate classical solutions which are more or less 
static in time besides non-static ones. 

For both kinds of boundary conditions, among the solutions we have found there 
are such which can be easily identified as CAL and DD. In order to allow a simple 
distinction between the non-static calorons CAL and the dissociated but static DD pairs 
we introduce a quantity which represents for an Euclidean configuration the mobility of 
the action density. For brevity, we call it non-stationarity 

S t = ^2^2\s(x,t;fM,u) - s(x,t- l;(i,u)\/S. (21) 

t,X fl<V 

The action density per plane s(x, £; /x, v) and the normalization factor, the total action S, 
are defined in fll~5|) . We have monitored how frequently objects with given S t are found at 
the lowest action plateaus. The histograms of S t look similar for both types of boundary 
conditions. For (3 = 2.2 they have a peak at St = 0.02 — 0.04 and a long tail for large St. 

We have convinced ourselves that a cut S t < 0.17 well separates -D-D-objects which 
are static with two well-separated maxima of the densities of the topological charge qt{x) 
and action <j(x)). For St > 0.17 the objects can be easily interpreted as CAL which are 
non-static, with an approximately 0(4) rotationally symmetric action distribution, with 
a single maximum of qt{x) and q(x) ) in 3D space. Both DD and CAL are showing a 
pair of opposite-sign peaks of the Polyakov loop. 

4.1 DD pairs 

For a special DD solution found with p.b.c, we show in Figs. [I] (a) and (b) two- 
dimensional cuts of the topological charge density qt(x) (a) and of the Polyakov loop 
distribution L(x) (b). The DD solution was obtained from an equilibrium configuration 
representing (3 = 2.2, i. e. the confinement phase. We clearly see the opposite-sign peaks 
of the Polyakov loop variable correlated with the equal-sign maxima of the topological 
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charge density. The boundary values of the Polyakov loop are slightly varying because 
they are not fixed here to a well-defined value. This is the only difference observed be- 
tween the two types of boundary conditions. In principle, for p.b.c. the holonomy could 
be arbitrary. What really happens to the asymptotic holonomy we will discuss in Section 
||. As a consequence the ratio of the action carried by the well-separated dyon constituents 
can take any value. 

For the same DD solution, Figs. p] (c) and (d) show the scatter plot of the 70 lowest 
complex Wilson fermion eigenvalues (|20"D for k = 0.14, both for time-periodic (c) and 
time-antiperiodic (d) boundary conditions for the fermion fields. In both cases we find 
one isolated low-lying real eigenvalue which can be related to a zero-mode of the zero-mass 
continuum Dirac operator. The corresponding (projected) eigenmode densities ^i(j(x) are 
drawn below, in Figs. |T] (e) and (f). They show a localization behavior as analytically 
proposed in Eqs. flU]) . For the time-antiperiodic b.c. the eigenmode is localized at the 
dyon exhibiting the negative peak of the Polyakov loop related to Taubes winding |55|] . 

For the given solution created on the lattice we have carried out a fit with the analytic 
formula |28| to reproduce the action density (0) . This has provided the parameter values 
x\ = (8,5,11), x 2 = (5,8,5) and to = 0.202. Fig. |2] shows one-dimensional projections 
of the same gauge field configuration together with the analytical results produced with 
the given fit parameters for the topological density, the Polyakov loop and the fermionic 
mode density with time-antiperiodic and -periodic boundary conditions according to the 
expressions (|T4]). The last two curves are predictions, rather than fits. There is an 
impressive agreement with the numerical shape of the fermionic zero-mode density. 

Gauge fixing to MAG we can search for the Abelian monopole content of the field 
configurations under inspection. We are interested in the positions of the world lines of 
monopole-antimonopole pairs. For static DD solutions we always find a pair of static 
(anti)monopoles with world lines coinciding with the centers of the dyons. All these fea- 
tures have been observed for DD-solutions irrespective of the spatial boundary conditions 
imposed in the process of cooling. 

4.2 CAL configurations 

In Fig. we show a typical CAL solution, with an approximately 4D rotationally invariant 
action distribution, obtained at (3 = 2.2 from cooling with periodic boundary conditions. 
The configuration possesses a large value of the non-stationarity 5 t . Again we plot 2D 
cuts for the topological charge density, for the Polyakov loop and the fermionic eigenvalues 
together with the eigenmode density for the distinct real eigenvalue. The full topological 
charge Q t is unity. The expected pair of narrow opposite-sign peaks of the Polyakov loop 
is nicely visible. 

The fermionic zero-modes for time-periodic and time-antiperiodic b.c. for this con- 
figuration are only slightly shifted relative to each other. A reasonable fit with the an- 
alytic solution can be obtained showing that this caloron is nothing but a limiting case 
of the generic DD solutions. At that point we may conclude that cooling, even with 
non-fixed holonomy, yields almost-classical solutions which show all characteristica of the 
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KvB calorons. The typical CAL configurations show, after putting them into MAG, a 
closed Abelian monopole loop circulating around the maximum of the action density in 
the 4Z}-space. 



4.3 DD pairs 



As previously observed for the case with f.h.b.c. [p4|| , we have found also other field 
configurations with an action on the instanton level, S ~ Si n st> which are very stable 
against cooling which motivates us to interprete them also as approximate solutions of 
the lattice equations of motion. With very low non-stationarity St = 0.004 ± 0.002, they 
consist of two lumps of action with opposite-sign topological charge densities. We call 
them dyon-antidyon pairs, DD. Each of their lumps turns out approximately (anti) self- 
dual. The total topological charge of the entire configuration is always zero. Therefore, 
each of the lumps carries half-integer topological charge. The Polyakov loop exhibits two 
peaks, in this case of equal sign. 

Also for DD pairs MAG fixing offers an Abelian monopole interpretation. After 
Abelian projection a static Abelian monopole-antimonopole pair can be found at the 
positions of the topological charge centers. 

Searching for the eigenvalues of the Wilson-Dirac operator we did not find real eigen- 
values but sometimes pairs of complex conjugated eigenvalues with very small imaginary 
parts. This feature is very similar to that of dilute superpositions of instantons with 
antiinstantons. This gives us good reason to interpret DD pairs as superpositions of 
single BPS solutions with half-integer topological charge. To the best of our knowledge, 
analytic solutions of this kind have not been reported in the literature. An example for a 
DD pair is reproduced in Fig. f|. We did not find any real or near-to-real modes for the 
time-antiperiodic boundary conditions. 



4.4 The composition of the lowest action plateaus 

Finally, by cooling with both kinds of spatial boundary conditions we have found objects 
becoming very stable at even lower action, for lattice size 16 3 x4, S ~ Sinst/2 and ~ Si nst /4. 
Their (color-) electric contribution to the action is very small compared with the magnetic 
contribution. Moreover, they are perfectly static with St = 0.003±0.002. Employing MAG 
we have convinced ourselves that they are purely Abelian. In the confinement phase they 
are occuring quite rarely directly in the cooling process. They are more common to appear 
after monopole-antimonopole pairs observed at S ~ S^t annihilate in the final stage of 
the relaxation. Therefore, we shall not consider them in detail here. But they seem to play 
an important role in the deconfinement phase [|37[ |35| . Since they are purely magnetic 
solutions - pure magnetic fluxes or 't Hooft-Polyakov-like monopoles - we shall abbreviate 
them M. 

In Table 1 the relative frequencies to find different types of classical configurations 
(DD, CAL, DD and M) at and below the one-instanton action plateaus are shown. We 
compare here p.b.c. with f.h.b.c. For each j3 < (3 C we have investigated O(200) Monte 
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Carlo equilibrium configurations. We can conclude that cooling applied to configura- 

Table 1: Relative frequencies of the occurrence of different kinds of (approximate) solu- 
tions, for different values of (3 and for both kinds of boundary conditions of the gauge 
field (f.h.b.c. and standard p.b.c). The lattice size is 16 3 x 4. 



Type of 
solution 


Boundary 
condition 


P = 2.20 


P = 2.25 


P = 2.30 ~ p c 


DD 


f.h.b.c. 
p.b.c. 


0.46 ±0.05 
0.43 ±0.05 


0.52 ±0.05 
0.44 ±0.05 


0.45 ±0.05 
0.23 ±0.03 


CAL 


f.h.b.c. 
p.b.c. 


0.19 ±0.03 
0.24 ±0.03 


0.17 ±0.03 
0.26 ±0.03 


0.15 ±0.03 
0.26 ±0.03 


DD 


f.p.b.c. 
p.b.c. 


0.28 ±0.04 
0.18 ±0.03 


0.26 ±0.04 
0.16 ±0.03 


0.31 ±0.04 
0.10 ±0.02 


M 


f.p.b.c. 
p.b.c. 


0.01 ±0.01 
0.04 ±0.02 


0.01 ±0.01 
0.03 ±0.01 


0.03 ±0.01 
0.10 ±0.02 


trivial vacuum 


f.p.b.c. 
p.b.c. 


0.06 ±0.02 
0.11 ±0.02 


0.04 ±0.02 
0.11 ±0.02 


0.06 ±0.02 
0.31 ±0.04 



tions in the confinement phase produces all objects with relative probabilities which are 
approximately independent of the type of boundary conditions imposed. 

For the deconfinement phase we have seen that the strong enhancement of DD config- 
urations earlier found for cooling with f.h.b.c. [34| (which would be compatible with the 
suppression of the topological susceptibility) is not reproduced under cooling with stan- 
dard p.b.c. In the standard case, the probability to obtain any topologically non-trivial 
object drops sharply with P > P c . Cooling down to the one-instanton action plateaus pro- 
vides only trivial vacuum or M-configurations. Because this latter observation was based 
on a physically small 3-volume (16 3 *4 for P = 2.4) finite-size effects might have been too 
strong to preclude a final conclusion. The structure of cooled deconfined configurations 
will be addressed in a further investigation. 

The independence of the boundary conditions, however, in the confinement phase has 
to be taken seriously: the enforcement of a L = boundary condition seems to be not far 
from the situation with standard periodic boundary conditions in the MC equilibrium. 
Some details will be discussed in the next Section. 



5 Dilute gas configurations of dyons and antidyons 
at higher action plateaus 

Within the confinement phase, for < T < T c and for both kinds of spatial boundary 
conditions, we have also studied in detail semi-classical configurations at higher action 
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plateaus. They represent snapshots of earlier stages of the cooling histories because the 
stopping criteria were focussed on multiples of the instanton action. This study should 
allow us to observe superpositions of classical solutions studied in Section ^ promising to 
be relevant for a semi-classical approximation of the non-zero T partition function. So far 
in the literature, the semi-classical approach to QCD at non-zero temperature is entirely 



based on Harrington-Shepard caloron solutions with trivial holonomy |27], 26]. Our main 
concern is here, whether superpositions of solutions with non-trivial holonomy naturally 
occur under cooling. 



5.1 Landscapes of topological density and Polyakov loop, fermion 
zero modes and Abelian monopoles as dyon finder 

Therefore, we expose equilibrium Monte Carlo lattice gauge field configurations to cooling, 
this time stopping under criteria which apply to different, subsequent action windows. We 
have been monitoring the landscape of topological density, of the Polyakov line operator 
as well as the localization of the fermionic zero-modes in the semiclassical candidate 
configurations. 

Searching for more complex approximate classical solutions we have modified our 
previous stopping criterion, triggering now on: 

• {m-\)S imt <S n < (m+±)S- mst , m = 2,3, and 

• S n — 2 S n -i + S n -2 < 0. 

In particular we inspected the first (highest) visible plateaus which occurred at various 
m-values, typically in the range m ~ 8 4- 20 for (3 = 2.2 on a lattice of the size 16 3 x 4. 
Then, we looked at the series of subsequent action plateaus. In terms of the objects 
classified in Section f|, we have scanned the resulting plateau configurations. 

For more definiteness concerning the moment of taking snapshots of the configurations 
undergoing cooling along a plateau, we have additionally introduced a measure A for the 
mean violation of the classical lattice field equations per link (see ||32|| ) 

A = ^ v ^E{ tr ^-^)(^-^) t ]} 1 , (22) 

where 

U x ,fi = c y ] ^Ux^Ux+u^U^^ j; + Ul^z^Ux-v^Ux+ft-u^u (23) 

is the local link x, fj, being the solution of the lattice equation of motion, with all degrees 
of freedom coupled to it being fixed. The factor c is just a normalization of the staple 
sum such that U Xtfl G SU (2) Q 

On the first visible plateau we find a gas of localized lumps carrying topological charge, 
where an identification in terms of dyons D and/or antidyons D is still difficult. 

1 The replacement U x ,p, — ¥ U x .^ is exactly the local cooling step as applied throughout this paper. 
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Independent of the kind of boundary conditions employed, at somewhat lower action 
plateaus with m < 10, we are able to clearly recognize dyons D and antidyons D carrying 
non-integer topological charges. During cooling more and more of these objects disappear. 
However, at all plateaus we observe an even number of peaks of the spatial Polyakov loop 
landscape For illustration see Figs. || - which show one and the same gauge 

field configuration at different stages of the cooling process. In this case f.h.b.c. have 
been used. 

In Fig. |5] we show the action (in units of the instanton or caloron action), the non- 
stationarity St, and the mean violation A of the lattice field equations per link. At three 
subsequent, already lower action plateaus (labelled by m) we indicate the iteration steps 
A (for m = 4), B (for m — 3) and C (for m — 2), respectively, where A passes through 
local minima. The corresponding semi-classical field configurations are then displayed in 
Fig. ^| by means of the 2.0-projected (by summing with respect to the third coordinate) 
spatial topological charge density and the 2D projected Polyakov loop distribution. More 
or less well one can recognize in these Figures that at stage A we have a superposition of 6 
dyons and 2 antidyons. The topological charge sector has been independently determined 
to be Q t = 2. A DD pair decays or annihilates from A to B such that we have 5 dyons 
and 1 antidyon at the next stage. The topological charge did not change. Finally, stage 
C exhibits a superposition of 4 dyons, again with Q t = 2. The latter configuration is 
very stable. While it stays at the same action over thousands of cooling sweeps, the 
non-stationarity S t gradually rises. A closer look then shows that the scale size of one of 
the dyon pairs shrinks, transforming this pair into a small undissociated and non-static 
caloron, which finally disappears after having turned into a tiny dislocation strongly 
violating the equations of motion (compare with Fig. |^). The average violation of the 
equations of motion per link peaks immediately before the configuration drops to the next 
plateau. 

This example shows that we have superpositions of non-integer Q t lumps, which can 
be interpreted as described in the previous Section. To make sure that this is really the 
case we provide also the eigenvalue scatter plots for the Wilson-Dirac operator for stage 
A (Fig. |7|), for stage B (Fig. |8|) and for stage C (Fig. ||). Fig. |7| shows four real fermion 
modes (under time-periodic boundary conditions) which could suggest an interpretation 
of configuration A as a superposition of 3 DD pairs and a single DD pair. However, the 
inspection of the time-antiperiodic b.c. case provides only two real modes which supports 
a dyonic content consisting of 2 DD pairs plus 2 DD pairs, an interpretation which 
naively is possible as well. In the stage C, also for time-periodic boundary conditions, 
we see clearly two real modes sitting on top of two dyon lumps. We have checked that 
the modes jump to the remaining dyon lumps when changing the fermionic boundary 
condition to time-antiperiodic. 

We have studied also the Abelian (anti)monopole structure after fixing to the MAG 
and Abelian projection. We see a strong correlation of the peaks of the Polyakov loop 
with the positions of the (anti)monopoles. This can be seen in Fig. |TD|. The pair structure 
in terms of Abelian monopoles, occuring on all action plateaus and the annihilation of 
single (monopole-antimonopole) pairs during further relaxation, provides an additional 
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signal for the topological content as superpositions of non-trivial holonomy CAL, DD 
(DD) or DD pairs. 



5.2 The role of the boundary conditions: are non-trivial bound- 
ary conditions "natural" ? 

In order to understand this from the point of view of single caloron solutions with non- 
trivial holonomy we have to find out whether (approximately) those asymptotic holonomy 
boundary conditions as typical for the dyonic (antidyonic) semi-classical background exci- 
tations are actually present during the cooling process when periodic boundary conditions 
(no particular holonomy boundary conditions !) are applied to the full volume. Then it 
would be easier to accept that similar (albeit fluctuating) boundary conditions might also 
be realized in the semi-classical vacuum. 

To answer this question we define the asymptotic holonomy of a cooled configu- 
ration as the average of L(x) over all points x in 3D space for which the 3D projected 
action density satisfies q(x) < .0001, i.e. it takes minimal values which are typically seen 
in deep valleys around the topological lumps. 

In Fig. |TT] we show histograms of obtained at different plateaus during the cooling 
histories of an ensemble of 0(200) configurations produced at (3 — 2.2 on a 16 3 x 4 
lattice with standard p.b.c. We see a clear peak at = for higher action plateaus. 
The distribution is more narrow, than the pure Haar measure would tell us. However, 
approaching lower- lying plateaus, the real distribution becomes flat. 

Closing the discussion of the local boundary conditions, let us finally concentrate on 
those configurations which belong to the bins of —1/6 < < +1/6, i.e. those which 
realize « (L) in the confined phase. For these configurations we turn our attention 
to the local correlation between the Polyakov loop L(x) and the action density values 
q(x) measured at the same positions x. We studied conditional distributions P[L\q] of the 
Polyakov loop values L(x) at spatial points where the spatial action density q(x) equals q. 
The distributions are normalized for each bin of q. The distributions are shown as surface 



plots in Fig. |12| for different plateaus. They show that with higher action density the 
corresponding local Polyakov loop values tend to be closer to the peak values L — ±1. 
The distributions do not depend on which action plateau they were collected. We have 
identified in Section |] KvB solutions and DD configurations at the one-instanton action 
plateau, and we have strong indications that the same objects occur on the higher action 
plateaus, too, i.e. in superpositions of topologically non-trivial lumps of action. This 
argument is also supported by a comparison with the same kind of conditional distribution 
obtained from the analytic KvB solution with random parameter distribution which we 
have discretized on the lattice. The resulting distribution P[L|s] is shown in Fig. [13] (a) 
and compared with the distribution for calorons with trivial holonomy (b). 
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5.3 The Taubes rotation in many-dyon configurations 

Finally, it is interesting to analyse the relative orientation of the dyons in color space. In 
Section ||| we have described how the Taubes winding in a DD caloron could be detected 
in a gauge independent way. The analysis is also here, for real cooled configurations, 
simplified by the observation that at the center of a dyon both electric Ef and magnetic 
B? field strengths (which satisfy (anti) self-duality Ef = ±5? ) form orthogonal matrices 
in color (a) and space (z) indices. Thus, it suffices to consider three (i = 1, 2, 3) vectors in 
color space Ef(x) forming a local reference frame at the center of a dyon with which any 
other local reference frame can be compared. The comparison can be made in a gauge 
invariant manner by connecting the centers at which the field strengths are measured by 
a parallel transporter, the Schwinger line. 

The exploration of a few of semi-classical lattice configurations containing superpo- 
sitions of several D and D has shown the following common features. All dyons with 
negative peak value of the Polyakov loop have more or less random color orientation rel- 
ative to each other, but this relative orientation is static along the time axis. For dyons 
with a positive peak value of the Polyakov loop holds the same. Also their orientation 
relative to each other seems to be random but static. But the relative orientation be- 
tween the L = +1 and the L = — 1 dyons is changing along the Euclidean time and the 
change is nothing but a homogeneous rotation in color space with the angle period 2tt. In 
the analytic solution representing just two dyons (or two antidyons) this rotation is also 
present. It has been discussed and related to Taubes winding in Sect. ||. 

To illustrate this observation let us consider a configuration obtained on the lattice 
of size 16 3 x 4 with f.h.b.c. after 400 cooling steps. The configuration contains three D 
and one D as shown in Table 2. We have taken the first D as the reference point with 

Table 2: Configuration of 3 D and 1 D, obtained after 400 cooling steps with f.h.b.c. from 
an equilibrium configuration produced at (3 = 2.2 on a 16 3 x 4 lattice. 



Type 


3D position 


Polyakov loop 


D 


x x = ( 5, 3, 2) 


-1 


D 


x 2 = (12, 6, 3) 


-1 


D 


f 3 = (13, 14, 5) 


+ 1 


D 


£4 = ( 5, 7, 13) 


-1 



respect to which the relative color orientations of the other dyons and the antidyon were 
determined. We have calculated the matrices (cf. eq. fll3p) Rl%(t), R\l{t) and R}^{t). 
They show the relative orientations of the objects n = 2, 3, 4 - represented by their electric 
fields E k (x n ,t) = E^Xn) 2 ^- which appear parallel transported to the position of the first 
object in the form S(xi,t; x n , t) E k (x n , t) S + (xi,t; x n , t) - with respect to the first object 
- represented by its electric fields E k {xi). Then the evolution in time of the relative 
orientation can be investigated. While the matrices R\lit) = tr (Ei(xx, t)SEk(x2, t)S + ), 
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and Rjk(t) turned out to be constant in time, the matrix R\lit) performed a color rotation 
with a constant angle increment n/2 from one time slice to the next time slice about the 
color axis n = (—0.919,0.278,0.275). The orientation of this axis seems to be random, 
but the rotation angle is well-defined. For N t — 4 it corresponds to a full color rotation 
over the full Euclidean time period. Our general observation illustrated by this example of 
a moderately complicated superposition shows that also these more complicated objects 
exhibit a strong correlation in the color orientation analogous to that present in a single 
DD KvB pair. A semi-classical approximation of the path integral would have to take 
into account this kind of color correlation. 

6 Conclusions 

We have generated SU(2) lattice gauge fields at non-zero temperature in the confinement 
phase. We have cooled them in order to analyse their topological content. Fixed holon- 
omy spatial boundary conditions have been used as well as standard periodic boundary 
conditions. The results for these two kinds of boundary conditions semi-quantitatively 
agree with each other. This is specific for ensembles describing the confinement phase. 

Independently of the boundary conditions we have found superpositions of calorons, 
dyons and antidyons, the latter with positive and negative non-integer topological charges. 
The topological lumps appear also as peaks of the Polyakov loop modulus |L(x)|, with 
calorons being a limiting case with a close pair of L(x) = ±1. Investigating also the 
localization behavior of the real eigenvalue modes of the (non-Hermitean) Wilson-Dirac 
operator we could present convincing evidence that for calorons and for dyon-dyon pairs 
an interpretation in terms of KvB solutions is natural. Chosing time-periodic and time- 
antiperiodic boundary conditions for the fermions focusses on L(x) = — 1 or L(x) = 1 
dyons, respectively. 

On higher action plateaus we have found that the dynamics generically leads to non- 
trivial holonomy outside the lumps of action and topological charge. The multi-dyon- 
antidyon structure can be finally resolved by a combination of two tools: localization of 
fermionic real modes and the relative color-orientation of the color-electric field strength. 

A semi-classical treatment of the path integral at non-zero temperature close to the 
deconfinement phase transition should be built on superpositions of calorons, dyons and 
antidyons, however with the holonomy as a free parameter. 

We have seen that such superpositions would imply a strong correlation in the relative 
color orientation (and its Euclidean time dependence) between pairs of seemingly inde- 
pendent topological lumps. To the best of our knowledge, such superpositions have not 
yet been constructed analytically. 

It is already clear that the development of a semi-classical approach based on solutions 
with non-trivial holonomy is much more complicated than the instanton (caloron) gas or 
liquid, and it might turn out not to lead to a practicable model. 

Nevertheless, facing the non-trivial structures found in this paper might contribute 
to a better understanding of the mechanism driving the deconfinement transition. In as 
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far a certain working picture of a dilute gas of these configurations can be developed and 
whether it will improve our understanding of quark confinement itself remains an open 
question to which we hope to come back in the next future. 
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Figure 1: Various portraits of a selfdual DD pair obtained by cooling under periodic 
gluonic boundary conditions. The sub-panels show: appropriate 2D cuts of the topological 
charge density (a) and of the Polyakov loop (b), the plot of lowest fermionic eigenvalues 
(c,d) and the 2D cut of the real-mode fermion densities (e,f), for the cases of time-periodic 
(c,e) and time-antiperiodic (d,f) fermionic boundary conditions, respectively {(5 = 2.2 and 
lattice size 16 3 x 4). 
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Figure 2: The DD lattice configuration of Fig. f (dashed lines), fitted by the KvB 
solution (solid lines) according to the action density, shown in three spatial views (ID 
projections). In (af ),(a2),(a3) the topological charge density is summed over the two 
unspecified spatial coordinates. In (bl),(b2),(b3) the Polyakov loop is averaged over the 
two unspecified spatial coordinates. In the time-antiperiodic case (cl,c2,c3) or the time- 
periodic case (dl,d2,d3), respectively, the fermion density it is summed over the two 
unspecified spatial coordinates. 
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Figure 3: Various portraits of a self dual CAL configuration obtained by cooling under 
periodic gluonic boundary conditions. The sub-panels show: appropriate 2D cuts of the 
topological charge density (a) and of the Polyakov loop (b), the plot of lowest fermionic 
eigenvalues (c,d) and the 2D cut of the real-mode fermion densities (e,f), for the cases of 
time-periodic (c,e) and t ime- ant iper iodic (d,f ) fermionic boundary conditions, respectively 
(/3 = 2.2 and lattice size 16 3 x 4). 
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Figure 4: Various portraits of a mixed-duality DD pair obtained by cooling under 
periodic gluonic boundary conditions. The sub-panels show: appropriate 2D cuts of the 
topological charge density (a) and of the Polyakov loop (b), the plot of lowest fermionic 
eigenvalues (c,d) for the cases of time-periodic (c) and time-antiperiodic (d) fermionic 
boundary conditions, respectively ((3 = 2.2 and lattice size 16 3 x 4). A 2D cut of the 
fermionic mode density related to the two distinct almost real eigenvalues in (c) is shown 
in (e). 
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Figure 5: Part of the cooling history for a gauge field configuration taken from the Monte 
Carlo sample generated at (3 = 2.2 on a 24 3 x 4 lattice, with f.h.b.c. of L £en = 0. The 
sub-panels show: (a) full action S/Si nst , (b) non-stationarity S t and (c) mean violation A 
per link of the lattice field equations, vs. the number of cooling steps. The vertical dotted 
lines indicate the passages of A through local minima having occured at 800 (A), 1650 
(B) and 7000 (C) cooling steps for which the configurations will be portrayed in Figures 
§, and g 
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Figure 6: Configurations on the 24 3 x 4 lattice (from equilibrium at j3 = 2.2 with f.h.b.c.), 
as indicated in Fig. § after 800 (A,A'), 1650 (B,B') and 7000 (C,C) cooling steps. (A,B,C) 
show 2D projections of the topological charge density qt(x) and (A',B',C) of the Polyakov 
loop L(x), respectively. Cooling has been employed with f.h.b.c, Lg £n = 0. 
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Figure 7: The lattice field configuration depicted in Fig. || (A, A') for 800 cooling steps 
(f.h.b.c). Here (A) plots the eigenvalues of the Wilson-Dirac operator in the complex 
plane for k = 0.140 and the case of time-periodic fermionic b.c; (Al,...,A4) show 2D 
projections of the fermionic mode densities related to the four distinct real eigenvalues. 
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Figure 8: The lattice field configuration depicted in Fig. ^| (B,B') for 1650 cooling steps 
(f.h.b.c). Here (B) plots the eigenvalues of the Wilson-Dirac operator in the complex plane 
for k = 0.140 and the case of time-periodic fermionic b.c; (B1,B2) show 2D projections of 
the fermionic mode densities related to the two distinct almost real eigenvalues, whereas 
(B3,B4) present the densities related to the two real eigenvalues. 
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Figure 9: The lattice field configuration depicted in Fig. ^ (C,C) for 7000 cooling steps 
(f.h.b.c). Here (C) plots the eigenvalues of the Wilson-Dirac operator in the complex 
plane for k = 0.140 and the case of time-periodic fermionic b.c; (CI, C2) show 2D 
projections of the fermionic mode densities related to the two real eigenvalues. 
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Figure 10: Histogram P(L) of the values of the Polyakov loop L(x) taken at x where 
time-like Abelian (anti)monopoles are found. The data represent cooling plateaus at 
m = 4 obtained at f3 — 2.2 with lattice size 16 3 x 4. 0(2400) non-vanishing monopole 
charges were collected. 
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Figure 11: Histogram P(L 00 ) of the values of the Polyakov loop at "infinity" (as explained 
in the text) seen on the first plateau (a) and at plateaus with m ~ 4 (b), m ~ 3 (c), 
m ~ 2 (d), m ~ 1 (e). For comparison the distribution expected from the pure Haar 
measure Puaa,r(L) ~ Vl — L 2 is shown with the same normalization (dashed lines). The 
equilibrium ensemble was generated at f3 — 2.2, the lattice size is 16 3 x 4, cooling was 
performed using periodic b.c, 0(200) configurations were investigated. 
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Figure 12: Conditional distributions P [L\q] relating local values of L(x) with the spatial 
action density q(x) for cooled configurations at plateaus with 7 > m > 5 (a), m ~ 4 (b), 
m ~ 3 (c) and m ~ 2 (d). The equilibrium ensemble was generated at /3 = 2.2, the lattice 
size is 16 3 x 4, cooling was performed using p.b.c. 
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Figure 13: Conditional distributions -P[X|<r] as in Fig. [T2L obtained for random KvB 
solutions (DD or CAL) (a) and for calorons with trivial holonomy (b), for comparison. 
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